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. A stochastic approach to the quantum dynamics randomly modulated in time by a discrete state 

' non-Markovian noise, which possesses an arbitrary non-exponential distribution of the residence 

times, is developed. The formally exact expression for the Laplace-transformed quantum propa- 
gator averaged over the stationary realizations of such N-state non-Markovian noise is obtained. 
' The theory possesses a wide range of applications. It includes some previous Markovian and non- 

' Markovian theories as particular cases. In the context of stochastic theory of spectral line shape and 

, relaxation, the developed approach presents a non-Markovian generalization of the Kubo- Anderson 

. theory of sudden modulation. In particular, the exact analytical expression is derived for the spectral 

' line shape of optical transitions described by a Kubo-oscillator with randomly modulated frequency 

which undergoes jump-like non-Markovian fluctuations in time. 
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'j^ ; I. INTRODUCTION 

> \ 

^ , The dynamics of classical and quantum systems in the presence of randomly fluctuating micro-fields of environment 
presents one of the fundamental problems in physics. Spin relaxation in solids (IJ, |3] , exciton transport in molecular 
^ \ systems , single- molecular spectroscopy , classical transport processes with fluctuating barriers present a few 
relevant examples. A popular approach consists in modeling the ambient noise influence by means of a classical 
stochastic field acting on the dynamical system. In the case of quantum systems, such a phenomenological approach 
is known under the label of Stochastic Liouville Equation (SLE) approach 1^ is suitable in the limit of 

sufficiently high (formally infinite) temperatures |^ . In the field of chemical kinetics a similar methodology is 
known under the label of rate processes with dynamical disorder where the rates of chemical reactions fiuctuate Q . 
' Moreover, the addition of a non-equilibrium classical noise into dissipative quantum dynamics can serve to describe 
^ \ the influence of the non-equilibrium environmental degrees of freedom on the transport properties [l^ . 
\^ . The tractability of the SLE approach, which allows to arrive at the exact model solutions in several known cases, 
^\ ' makes it popular over the years. For example, the case of two-state quantum dynamics subjected to a white Gaussian 
lO . noise can be treated exactly and the corresponding exact master equations for the averaged parameters of quantum 
' systems can be derived However, already in the case of a colored Gaussian noise (the so-called colored noise 

' problem) a perturbation theory must be used which leads generally to an approximate description, e.g., within a 
generalized master equation approach [T]|. 
— ^ ^ What to do, however, when the ambient noise has a non-Gaussian statistics and/or it does exhibit long-range, e.g., 
"j^ power law temporal correlations with a very large, practically infinite range? A relevant example is given by a 1//" 
' noise which is ubiquitous in the amorphous solids and other glass-like materials like proteins (12.. 13]. Any perturbation 

1 . theory in such situations will certainly fail and we are confronting with a rather difficult problem. Nevertheless, for 
'■^ ' an arbitrary quantum dynamics the problem of finding the corresponding noise-averaged propagator can be solved 

Ch , exactly, at least formally, for a rather general class of non-Markovian jump processes modeled as a continuous-time 
O ' random walk (CTRW) (14lll5lll6j within arbitrary, but finite number of states N . A similar problem has been already 
considered in several previous works, notably in Ref. [T^ llSj. These works did not solve, however, the problem at 
hand for the case of stationary noise-averaging for a multi-state non-Markovian noise, when the evolution of the 
considered stochastic process starts in the infinite past, and not simultaneously with the evolution of the considered 
dynamical system. In the case of a non-Markovian noise this problem is not trivial. Physically this is but the most 
' important and relevant case to study. 

In this work we utilize the most general description of the discrete state non-Markovian processes of the continuous 
time random walk (CTRW) type with uncorrelated jumps. Generally, such processes are defined by the set of 
probability densities ipij (r) for making transitions among the discrete noise states . The noise-averaged quantum 
propagator is obtained below for this general case. The approaches of Ref. ^3 ^''^d Ref. are unified within this 
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general description. Next, the problem of stationary noise averaging is considered. It is shown in a constructive way 
how to make the stationary noise averaging in the case of factorized probability densities like in Ref. |l8l |. but with 
the time-independent matrix of transition probabilities pij . In this respect, the present work is most close to Ref. jl9j , 
where a similar approach has been proposed, however the explicit solution for the stationary noise averaged quantum 
propagator has not been obtained. The corresponding formally exact expression for the Laplace-transformed averaged 
propagator is found in this work in the explicit form for the first time. This presents the first main result of this work 
which possesses an ample range of applications. In particular, the noise-averaged propagator of the Kubo-oscillator 
with a stochastically modulated frequency, as well as the corresponding line shape form are also obtained. This 
presents the second main result of this work which can be important, e.g., for the single-molecular spectroscopy. 



II. MODEL AND THEORY 

Let us consider an arbitrary quantum system with the Hamilton operator H[^{t)] which depends parametrically 
on the stochastic process ^(t) which in turn can acquire randomly in time N discrete values ^i. Accordingly, the 
Hamiltonian H[£^{t)] can take on iV different random operator values H[^i]- The discrete stochastic process is assumed 
to be a non-Markovian renewal process which is fully characterized by the set of probability densities ipijir) which 
describe random transitions among the states ^i. Namely, ipijiT) is the probability density for making transition 
from the state j to the state i. These probability densities must obviously be positive and obey the normalization 
conditions 



N 

E 
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V-.j-Mdr^l, (1) 







for all j ~ 1, N. All random transitions are assumed to be mutually independent [2l|. The residence time distribution 
(RTD) ipjir) in the state j reads obviously 

i 

The survival probability (r) of the state j follows then as 

POO 

'^jir) = / ^,{r)dT. (3) 



This is the most general description used in the CTRW theory |l.6j. 

The problem is to average the quantum dynamics in the Liouville space which is characterized by the Liouville- 
von-Neumann equation 

= ~^Cm]pit) (4) 

for the density operator p{t) over the realizations of noise ^{t). C[^{t)] in Eq. Q stands for the quantum Liouville 
superoperator, C[£,{t)]{-) = j:[H[^{t)], (•)]. In other words, one has to find the noise-averaged propagator 

{U{t)) = (Texph* r^'mT)]dT]), (5) 

Jto 

where T denotes the time-ordering operator. It is the major advance of this work that we obtain the Laplace-transform 
of propagator ij^l in the exact form for arbitrary non-Markovian processes ^(t) of the discussed form. The results of 
previous research work done within the framework of SLE approach for the discrete Markovian processes j2j follows 
as a particular Markovian limit of the developed non-Markovian theory. Moreo ver, some previous non-Markovian 
theories, notably that by van Kampen and that by Chvosta and Reineker are also included as different 
interpretations of this general formulation. 

In particular, the approach of van Kampen is reproduced by introduction of the time-dependent age-specific rates 
kij (t) like in the renewal theory E^l • The probability densities ipij (r) then read 0, 

V'jj(T) fcy(T)exp[-'V' / kij{t)dt]. (6) 
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The Markovian case corresponds to kijir) = const. Any deviation oiipij{T) from the strictly exponential form which 
yields a time-dependence of the transition rates kijlr) amounts to a non-Markovian behavior j28 |. Furthermore, the 
survival probability $j (r) in the state j within the time-dependent rate description is given by 

N 

hj{t)dt] (7) 



^j{t) = cxp[-^ / 
and Eq. © then can be recast as 



V-yXr) :=%(t)$j(t). (8) 

The introduction of time-dependent rates is one possible way to describe the non-Markovian effects It is not 

unique. For example, Chvosta and Reineker adopted a quite different and more general standpoint Il8|. Namely, they 
defined 

i^tji-r) ■.^P^J{T)i^J{T) (9) 

with 'YliiPiji'^) = 1- The interpretation is as follows. The process stays in the state j during a random time interval 
characterized by the probability density ipj^r). At the end of this time interval the process makes jump into the state 
i with a generally time-dependent conditional probability pij (r) . Indeed, any stochastic process of the considered 
kind can be interpreted in this way. For some particular applications in [T^ the probability densities V'j (''■) were taken 
strictly exponential and all the non-Markovian effects were assumed to come from the time- dependent transition 
probabilities pij^r). By equating Eq. ((SJ and Eq. ^ and taking into account that ^/'j(r) := —d^j{T)/dT it is easy to 
see that van Kampen approach can be reduced to that of Chvosta and Reineker with the time-dependent transition 
probabilities 

and with the non-exponential probability densities ipjir) which follow as ipj{T) = exp[— JJ" ^j[t)dt\ with 7j(r) := 

The description of non-Markovian effects with the time-dependent transition probabilities pij (r) seems, however, be 
merely a theoretical device. It appears to be rather difficult (if possible) to obtain (r) from the sample trajectories 
of an experimentally observed random process f (t). In view of Eq. I|1U|) the same is valid for the concept of time- 
dependent rates. These rates cannot be measured directly from the sample trajectories. On the contrary, the RTD 
ipjir) and the time-independent pij (withp^i := what is assumed in the following) can be routinely deduced from the 
sample trajectories measured in a sm(7fe-mo/ecM/ar experiment. This latter description is definitely more advantageous 
from the practical point of view. From the experimentally well-defined quantities, ^jir) and Py , the corresponding 
time-dependent rates description can be readily found as 

, , , dln[$,(T)] , , 

hj{r)^-p., (11) 

Moreover, as it will be shown below namely this description with constant pij in Eq. Q does provide the consistent 
way in order to construct the stationary realizations of the non-Markovian process ^{t) and, therefore, in order to find 
{U{t)) averaged, correspondingly, over the stationary realizations of^{t). These circumstances give definite advantages 
of the approach with factorized ipijiT) ^^'^ time-independent pij as compare with the formulations in \1T^ and |l8j |. 
even though the technical details are quite similar. The correct quantum-mechanical propagator averaged over the 
stationary realizations of a non-Markovian discrete-state noise with an arbitrary number of states is obtained in this 
work for the first time. This corresponds to the situation where the quantum system was prepared at the time to in a 
non-equilibrium state described by the density matrix p(to) i but the noise was not specially prepared but it has been 
already relaxed to its stationary state. The ambient noise evolution was started in the infinite past and one assumes 
that the initial preparation of the quantum system in a non-equilibrium state has no influence on the noise source. 
This is the most important physical situation to confront with. 

The task of performing the noise-averaging of the quantum dynamics in Eq. jSJ can be solved exactly due the 
piecewise constant character of the noise ^ (t) [25l . Namely, let us consider the time-interval [to , t] and to take a 
frozen realization of ^{t) assuming k switching events within this time-interval at the time-instants ti. 



to < ti < t2 < ... <tk <t. 



(12) 
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Correspondingly, the noise takes on the values Cio'Cjn ■■■i^jk time sequel. Then, the propagator U(t,to) reads 

obviously 

U{t,to) = e"'''^['5^fcl(*"*'=)e"*^[^-*'=-il'*''-~*'=-i\..e~*'^[«^ol(*i~*») (13) 

Let us assume first that it is known with certainty that at the time instant to the process ^(t) has just started its 
sojourn in the state jo- In other words, the process ^(t) has been prepared in the state jo to- Then, the corresponding 
fc— times probability density for such noise realization is 

Pkiij, , tfc; , ifc-i; Cji , ti |Cio , ^o) = {t - tk)ijj,j,_, (tk - tk-i)...'4Jh]o (^i - *o)- (14) 

In order to obtain the noise-averaged {U(t,to))j„ conditioned on such nonstationary initial noise preparation one has 
to average (I13II with the probability measure in H14() (for k = 0, oo). Literally (operationally) this means the following. 
First, one has to construct the time-ordered product of p3|) and (|14|l . i.e. 

'l>,At~tk)e-^^^^^>'^^'-'''H,,,,_Ah (15) 

Second, one has to perform the fc-dimensional time integration of (|15|) over the variables {tk} within the time-ordered 
domain (|12|) and to sum the results over all possible {jk}- Furthermore, this procedure has to be repeated for every 
k = Q,N and the results summed at the end. The case fc = is special with 

Poiiio,h)^^,„{t-to). (16) 

The just outlined task can be easily done formally by use of the Laplace-transform, denoted in the following as 
A{s) := exp(— sT)A(T)dT for any time-dependent quantity A{t). For this goal, let us introduce two auxiliary 

matrix operators A[s) and B{s) with the matrix elements 

POO 

Aki{s):=5ki <^i{T)e-^'+'^^^'^'^^dT (17) 
Jo 

and 

POO 

Bki{s):^ Vfci(T)e-(^+^^[«'""dT, (18) 
Jo 

correspondingly. Then all what remains to do is to sum the geometric matrix operator series. The result reads, 

(C^(^))jo-E(^(-'^)[^-^Wr') : (19) 

where / is the unity matrix. It is quite obvious that instead of the quantum Liouville operator in Eq. Q there 
could be any linear operator, e.g. be a matrix and p(t) be then a vector function. Then, the developed theory 
can be immediately applied to the averaging of arbitrary linear stochastic differential equations |26l |2^ . Even 
nonlinear stochastic differential equations can be attempted to deal with by introducing a Liouville equation for the 
corresponding classical probability density ■ The earlier results in [l3| and are reproduced immediately from 
Eqs. (|17|I - H19|I by an appropriate modification of the considered problem and specifying the transition probability 
densities ■0^ (r) from a most general (non-factorized form) to a particular representation in accord with the above 
discussion. 

The derived result in Eas. H17|l - H19(l corresponds to the initial preparation of ^(t) in a particular state jo- Exper- 
imentally, this presents a quite unusual and strongly non-equilibrium situation. For a stationary environment one 
has to perform yet an additional averaging of {U{s))jg over the initial distribution pjg{to) taken as the stationary 
distribution, i.e., Pjoito) = where pf = liuit^oo Pjit)- Indeed, this presents a valid prescription for stationary 
noise-averaging in the Markovian case. However, in a non-Markovian case this prescription is not sufficient. 

Quite generally, the stationarity of noise realizations in the strict sense requires '27f that not only the single-time 
distribution Pj{t), but also all the multi-time joint probability distributions of the given process must be invariant 
of a simultaneous time shift of all time arguments. However, for many physical applications, the stationarity in a 
weak sense, i.e., on the level of the two-time, P{j,t; jo,to), joint distribution is sufficient. Then, e.g., the stationary 
power spectrum of the corresponding process can be defined. This two-time joint distribution can be expressed 
as P{j,t; jo,to) = Iijja{t\to)pjo{to) via the conditional probabilities Hjj(,(t|io) (propagator of the process). For a 
stationary process the consistency condition, p|' = X^jo Hjjo (^l^o)p|* , must be satisfied for all times, i.e. p** has 
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to be the fixed point of the corresponding propagator (see Appendix A). The propagator of non-Markovian process 
having this property can be called quasi-stationary. In the present case, in order to construct such propagator and the 
corresponding stationary realizations of the noise trajectories the probability density of the first time intervals must 
differ from the all subsequent ones. Indeed, if a noise state j was occupied at i = to with the stationary probability 
pj*, it is not known for how long this state was already occupied before to- The proper conditioning on and averaging 
over this unknown time must be made and the corresponding survival probabilities for the first residence time interval 
To = ti — to be introduced . These survival probabilities read [20l Isoj 

where (tj) = <^j{T)dT is the mean residence time (MRT) of the noise in the state j. Only for strictly exponential 
survival probabilities, i.e., in the Markovian case, '^^p {t) — <&j(T). Otherwise, this is not the case. The corresponding 

residence time distributions follow immediately as the negative time derivative of (r) in Eq. H2()|l and are known 
to be [1II23 

<V) = ^-i^ . (21) 

The first time transition densities ■(/',?'' (''') then follow as 



(22) 



For the logical consistency of this definition with the consideration pursued in remark |29l | Pij must he time-independent 
constants. Otherwise, a logical problem emerges: How to make the proper conditioning of pij{T) on the unknown 
times before to? This is the reason why within the approaches of time-dependent pij^r), or time-dependent rates 
kij (t) it is rather obscure how to solve the problem of stationary noise-averaging. Therefore, these approaches do not 
seem suit well for this stated purpose. 

In accord with the above discussion, the transition density ipjijo (^1 ~*o) in Eq. I|14|) must be replaced by V'jijo (^1 ~io) 
from Eq. (|^ . Moreover, ^jg(t — to) in Eq. lfTB|l must be replaced by $^-°^(t — to) from Eq. (^0)1 . To account for these 



3S A^^i^ (s) and (s) are introdi 



modifications, two auxiliary quantities A^j^j^ (s) and i?^^^ (s) are introduced which are given by the expressions similar 



to Eq. ((nil and Eq. but with $,^(t) instead of $j (t) and iplj^r) instead of i^tjir), i.e 

Ski 
in) 



<(^) ^ I <fKT)e-(^+^^[«'l)^dr = ^Ms). (24) 



and 

in) Jo in) 

The resulting geometric operator series can again be easily summed exactly. For the (stationary) noise-averaged 
Laplace-transformed propagator {U{s)) we obtain: 



{Uis))=J2[A<-'Hs)+A{s)[I-B{s)]-'m{s))^^pf, (25) 

where 

pf = (26) 

are the stationary occupation probabilities of (see Eq. (|A10|I in Appendix A and the corresponding discussion). This 
result can be brought into a physically more insightful form by using the identity <^j{T')dT' — (tj) — /J" ^j{T')dT' 
and upon introducing two new auxiliary quantities 

Cki{s):^Ski e-(^+^^[«'I)M <i>,(T')dT'dr (27) 
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and 

POO 

Dki{s):=5ki V'/(T)e-(^+^^[«'l)^dr . (28) 

Finally we obtain 

{U{s)) = {U{s)),,auc - ^ E - ^(^)[^ - PD{sr'PA{s))^^, (29) 

where {U {s)) static is the Laplace-transform of the statically averaged propagator 

(t/(r)).,at.c:=$]e-^^[«'^-l>r, (30) 

A: 

T = 'Yliki'^k) is the sum of mean residence times, and P is the matrix of transition probabilities pij ("scattering 
matrix" of the random process £,{t)). The result in Eqs. (|^ . (ISOJ together with Eqs. (fTTji . l(?7jl . presents the 

cornerstone result of this work which can be used in numerous applications. 

III. APPLICATION: KUBO-OSCILLATOR 

As a simplest practical example we consider the averaging of the so-called Kubo-oscillator 

x{t) = ioj[m]<t) . (31) 

This particular problem appears in the theory of optical line shapes, in the nuclear magnetic resonance and in 

the single molecular spectroscopy |^|. In Eq. w[C(t)] presents a stochastically modulated frequency of quantum 

transitions between the levels of a "two-state atom", or between the eigenstates of a spin 1/2 which are caused by 
the action of a resonant laser, or magnetic field, respectively. The spectral line shape is determined through the 
corresponding stochastically averaged propagator of Kubo-oscillator as 01 

^ - lim Re[U[iLO + e)] . (32) 

Note that the limit e ^ +0 in Eq. is necessary for the regularization of the corresponding integral in the 

quasi-static limit T oo. By identifying C[S,k] with —Uk in Eq. (|29|) we obtain after some algebra 



E 



1 - iIji{s - iuJi) ( I \ 1 - Ipnis - iiUn) 



with Dnm{s) ~ Snmi-'mis ~ i^m)- The Corresponding line shape follow immediately from Eq. (I33|l by virtue of Eq. 
()32|l . This result presents a non-Markovian generalization of the earlier result by Kubo 2] for arbitrary A^-state 
discrete Markovian processes. The generalization consists in allowing for arbitrary non-exponential RTDs ipkiT), or, 
equivalently, in accordance with Eq. Hll|) for time-dependent transition rates kij{T). This generalization is obtained 
here for the first time and presents one of our main results. Let us further simplify the result in Eq. H33() for the case 
of two-state non-Markovian noise with pi2 = P2i = 1. Then, Eq. ()33|l yields after some simplifications: 



{U{s)) = E 



{rk) 



k=l 2^~^^^ ^"^1^ + ^^2) 



{uJi-L02f [I - li)i{s - iuJi)\[l - Tp2{s - iuJ2)] ^^^^ 



(('^i) + {'^2)){s - iwiY{s - iW2Y 1 - ?/'i(s - iuji)%i)2{s - iuJ2) 



With H34|) in H32|l one obtains the result for the corresponding spectral line shape which is equivalent to one obtained 
recently in Ref. fsH using a different method. It is reproduced here as a simplest application of our more general 
approach. 
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IV. SUMMARY 

In this work the problem of the stochastic averaging of a quantum dynamics with non-Markovian fluctuating 
parameters has been investigated within the trajectory description of continuous time random walk theory. The 
formally exact expression for the stochastically averaged quantum-mechanical propagator is obtained for the most 
general CTRW with uncorrelated jumps and for a non-equilibrium noise preparation. The problem of stationary noise 
averaging has been solved for the practically relevant formulation with the time-independent matrix of transition 
probabilities pij . Especially, the formally exact expression for the stationary averaged quantum propagator has been 
found in an explicit form. This general expression has been used on order to find the stationary propagator of the 
Kubo oscillator describing the spectral line shape of optical transitions in a two-state atom. Further applications, 
such as decoherence of a two-state quantum dynamics driven by two-state non-Markovian noises, including 1//" noise 
case, are in progress. 
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APPENDIX A: PROPAGATORS AND GENERALIZED MASTER EQUATIONS 

In this Appendix, both the nonstationary propagator and the propagator for quasi-stationary initial preparations 
for the considered non-Markovian processes are obtained, along with the corresponding generalized master equations 
(GMEs). The propagator Hij of the process £,{t), or the matrix of conditional probabilities connects the initial 
probability vector p(to) with the final one, p{tQ + t), i.e., 

j 

The expression for the Laplace-transform 11^ (s) can be obtain in a way similar to the averaging of quantum propagator 
in Eqs. H19|) and (|25|l . Basically, one has to put there C ^ 0. For the nonstationary propagator of £,{t), i.e., when the 
process ^(t) starts its evolution at t = to in a, particular state the result reads for the general case of non-factorized 
probability densities i^ijij) follows: 

ri,,(s) = (A2) 

where 5* is the matrix of ijjij{s) and $j (s) = [1 — '!/'j(s)]/s, tpjis) = y^.- tl^i i (s). For quasi-stationary initial preparations 
the propagator of non-Markovian process is generally different [24. l32l|. In the factorized case ipijir) = pijipjij) it 
reads for the considered process (cf. Eq. (|25|l '): 

ng(.) -$r)(s)<5,, -I- pu,^, (A3) 

where $^°''(s) = 1/s — [1 — '0j(s)]/(s^(Tj)). The stationary populations follow as pf = Vm\s^Q{sIVij{s)). The quasi- 
stationary propagator n|j(r|0) must satisfy the consistency condition pf = nfj(T|0)p|* for all times r, i.e., pP* is 
the fixed point of nfj(r|0). Let us prove this fact and to find pf. It is more convenient to do both tasks by finding 
first the corresponding GMEs for Pi[t). These generalized master equations are of substantial interest per se. 

In order to find the corresponding GMEs the procedure of can be applied to a more general present case of 
non-factorized ipij{T)- Indeed, let us consider the conditional probability Pkij\ja)it) for making fc-jumps within the 
time interval [0,t] starting at < = in the state jo with the probability Pjo{0) and finishing in the state j with the 
probability Pj(t). This probability is given by a corresponding fc-dimensional integral of H14|l l (see the discussion 

below Eq. 11411 ) with the summation made over jk-i, jk-2i ■■■iji- The corresponding Laplace-transform Pjj^is) reads: 

PjSi') = E ••■E^^(^fe-i(^)--4..i(^)^S„(^) (A4) 

jk-l jl 
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for fc > 0. Here, Eq. H14() was used in a slightly modified form with ipjjH''') for the probability densities of the first time 

intervals. Furthermore, for fc = 0, Pjjg{s) — ^'f\s)Sjj„, where $j"^(s) is the Laplace-transform of the corresponding 

survival probability ^J^\t) = '^l^fj' {T')dT' of the first time-interval. For k > 2, The quantities Pjjg{s) satisfy 

obviously the following recurrence relation: 



pSi^)=^A^)j:^M^f^- (A5) 



Furthermore, the Laplace-transform of propagator (s) is expressed in terms of Pjj^ (s) as fljj„ (s) — J2h=o ^jja (*)■ 
Then, by virtue of Laplace-transformed Ea. (|Aip with to = 0, 



oo 



jo jo 



$.(^) E ^^"(^)r7T E E pS;'\')p^o{o) (A6) 

n ^"-y^) fc=2 30 



where the recurrence relation lA5p has been used. The use of X]fe°=2X]jo ^njo ^\^)P]oi^) — Pn{s) — $1"'' (s)p„(0) in 
Eq. (EEl) finally yields 

M^) = (^) E i44 + (^)p^- (0) 

^ - v;,„(s)^i^^ p„(0) . (A7) 

Let us consider now the case when the noise has been prepared at = in a particular state with the probability 
one. Then, tpf^\T) = ipijir) and the last term in IjATp vanishes. For this class of non-equilibrium initial preparations, 
the inversion of (|A7p yields: 



W(^) = - E / - t')pj{t')dt' f r,4t - t')pn{t')dt' 



(A8) 



where the Laplace-transformed kernels reads 



= (A9) 

with ipn{s) — ipjnis)- The just derived GME l)A8|l . (|A9|l is the most general GME for the continuous time random 
walk processes with uncorrelated jumps for the given class of initial preparations. In the case of factorized (but still 
non-separable) CTRW with ipij{T) = Pijipjir) it reduces to the GME of Burstein, Zharikov and Temkin |[19|. Moreover, 
in assumption that all ipjir) are equal (separable CTRW of MontroU and Weiss), it reduces further to the GME of 
Kenkre, MontroU and Shlesinger [33. The stationary populations can be obtained from (|A7|) as pf = \ims^Q{spj{s)). 

Assuming that the mean residence times (tj) exist, i.e. ipij{s) = — stij + o(s) with J^i'^ij = 1 ^^d ~ 
Eq. (jA7|) yields the system of linear algebraic equations for the stationary populations 



St ^st 

a.nv-^ (AlO) 



For an ergodic process the stationary probability to find the process in a particular state should be proportional to 
the time which the process spends in this particular state on average, i.e., be given by Eq. I|26(l . It is easy to verify 
that Eq. H2t)|l provides the solution of (|A10p for a,y — 1 (for the factorized case coincide obviously with pij). 
This latter condition can be also expressed as 



N 



l/'„(T)dT= 1 . (All) 
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For factorized case ipjniT) = Pjni'niT) with ')pj^\T) = ^j{T)/{Tj) |2^, the inversion of ljA7p yields the GME for the 
discussed non-Markovian process with quasi-stationary initial preparations: 

Pj{t) = -J Tj{t-t')[pj{t')-p^{0)]dt' +Y,PjnJ Tr^it ~ t')[p,,{t') - p^{0)]dt' 

with the kernels given by Tj{s) = sipj{s)/{l — tpj{s)). By choosing J3„(0) = p^* which satisfy Eq. IjAlOp it becomes 
obvious that Pn{t) — Pn is the solution of ljA12|l for all times t > 0. This means that the stationary jf* provides 
indeed the fixed point of the corresponding quasi-stationary propagator. 
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